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?^ ■ ABSTRACT 

We calculate the optical 6,, luminosity function of the 2dF Galaxy Redshift Survey 
Lj ■ (2dFGRS) for different subsets defined by their spectral properties. These spectrally 

r> I selected subsets are defined using a new parameter, 77, which is a linear combination of 

jrt ' the first two projections derived from a Principal Component Analysis. This parameter 

r] identifies the average emission and absorption line strength in the galaxy rest-frame 
spectrum and hence is a useful indicator of the present star formation. We use a 
total of 75,000 galaxies in our calculations, chosen from a sample of high signal-to- 
noise ratio, low redshift galaxies observed before January 2001. We find that there is a 
systematic steepening of the faint end slope (a) as one moves from passive (a = —0.54) 
to active (a — —1.50) star-forming galaxies, and that there is also a corresponding 
faintening of the rest-frame characteristic magnitude M* — 51og]^o(/i) (from —19.6 to 

— 19.2). We also show that the Schechter function provides a poor fit to the quiescent 
(Type 1) LF for very faint galaxies {Mtj — 51ogio(^) fainter than —16.0), perhaps 
suggesting the presence of a significant dwarf population. The luminosity functions 
presented here give a precise confirmation of the trends seen previously in a much 
smaller preliminary 2dFGRS sample, and in other surveys. We also present a new 
procedure for determining self-consistent X-corrections and investigate possible fibre- 
aperture biases. 

Key words: galaxies: distances and redshifts - galaxies: elliptical and lenticular, cD 

- galaxies: stellar content - galaxies: formation - galaxies: evolution 
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1 INTRODUCTION 

It is well established that the measurement of the galaxy 
luminosity function (LF) is sensitive to the type of galaxy 
being sampled. Morphologically early-type galaxies tend to 
be systematically brighter than their late-type counterparts, 
resulting in LF estimates with significantly brighter M* and 
shallower faint end slope a (e.g. Efstathiou, Ellis & Peter- 
son 1988; Loveday et al. 1992; Blanton et al. 2001). Similar 
trends are also present if one selects galaxies based on Ha 
equivalent widths (Loveday, Tresse & Maddox 1999), [Oil] 
equivalent widths (Ellis et al. 1996) or colour (Lin et al. 1996; 
Marzke & Da Costa 1997). Understanding and being able to 
quantify this variation in the LF is of great importance to a 
full understanding of galaxy formation and evolution. 

Up until now the greatest obstacle faced when attempt- 
ing to determine how the LF varied with different galaxy 
properties was the fact that one was required to divide al- 
ready small data sets, thereby losing much of the statistical 
significance in the LF estimations. In this paper we address 
this issue by making use of a subset of the galaxies observed 
to date in the 2dF Galaxy Redshift Survey (2dFGRS Mad- 
dox 1998; CoUess 2001). This subset is the largest single data 
set used in the calculation of LFs and can easily be divided 
several times whilst still maintaining very precise statistics. 

The 2dFGRS is a joint UK- Australian effort to map the 
distribution of galaxies down to an extinction corrected 6.7 
magnitude of 19.45 (median redshift 2 ~ 0.1). In so doing we 
expect to obtain 250,000 galaxy spectra from which redshifts 
can be determined. This is a factor of 10 more than any 
previous redshift survey. 

In order to divide our data set in a meaningful way we 
develop in this paper a classification of the galaxies based 
upon their observed spectra. As with other astronomical 
data, there are different approaches for analysing galaxy 
spectra (see e.g. Lahav 2000 and references therein). If one 
has a well defined physical model for galaxy spectra, then 
it is appropriate to estimate parameters of interest (such as 
age and star-formation rate) e.g. by Maximum Likelihood, 
directly using all the spectral bins, or via a compressed ver- 
sion of the data designed to give maximum information on 
the physical parameters of interest (e.g. Heavens, Jimenez & 
Lahav 2000). If, on the other hand, one prefers to let the data 
'speak for themselves', in a model-independent way, then it 
is more useful to look at the distribution of the galaxies in 
the high-dimensional space defined by the spectral bins. It 
is then possible to either look for distinct groups, e.g. early 
and late types, or more refined classes (e.g. Slonim et al. 
2000). An alternative, which we present here, is to find a 
continuous (sequence-like) parameterisation of the spectral 
features, which can later be divided into subsets. The pa- 
rameterisation we develop is denoted by 77 and essentially 
represents a measure of the average absorption or emission 
line strength present in each galaxy's spectrum. This clas- 
sification is robust to the known instrumental uncertainties 
and also has the advantage of being easily interpreted in 
terms of the current star-formation present in each galaxy. 

The set of galaxies we use in this paper comprises a 
sample of 75,000 galaxy spectra. This is more than 12 times 
larger than the one used in our previous analysis (Folkes et 
al. 1999, hereafter F99) and is by far the largest sample used 
to date. 



As well as showing the latest LF determinations we will 
also present a new procedure for calculating self-consistent 
A'-corrections and investigate possible fibre-aperture biases. 

The outline of this paper is as follows. Section H briefly 
summarises the 2dFGRS and describes the data set we are 
using. Section outlines our method for dividing the data set 
based on our parameterisation of the galaxy spectral type, 
r). Section gives a detailed description of the calculations 
involved in estimating the LF and in Section we discuss 
our results and future work. 



2 THE 2dFGRS DATA 

The 2dFGRS has already observed approximately 200,000 
unique galaxies for which it has obtained redshifts. The sur- 
vey, once complete, will cover approximately 20000° on the 
sky, split between two independent strips: one in the north- 
ern Galactic hemisphere and the other centred roughly on 
the south Galactic pole. In addition to this there are 99 
random flelds in the southern Galactic cap. 

The 2dF instrument itself is capable of observing up 
to 400 galaxy spectra simultaneously (Lewis et al., 2001). 
Each galaxy has been selected from a revised and extended 
version of the Automated Plate Measuring (APM) galaxy 
catalogue (Maddox et al. 1990) in order to determine its 
position and magnitude. It is then observed as part of the 
survey through a 140/im (~ 2 arcsec) diameter optical flbre. 
The fej magnitudes we use in this analysis are total mag- 
nitudes derived from updated versions of the original APM 
scans. These magnitudes have been updated to take into ac- 
count new GCD calibration data (see CoUess et al. 2001 and 
Norberg et al., 2001) and are believed to have an rms error 
of approximately ±0.15 magnitudes. 

Each of the spectra observed using the 2dF instrument 
spans 1024 channels with a spectral scale of 4.3A per pixel; 
the FWHM is measured by arc lines to be of the order of 1.8- 
2.5 pixels. Typically at the survey limit the observed spec- 
tra have an average signal-to-noise ratio of ~ 10 per pixel, 
sufBcient for determining redshifts and performing spectral 
analyses. For the purposes of our spectral analysis we have 
restricted ourselves to the redshift range of 0.01 < z < 0.2 
so that we are left with a uniform sample of galaxy spectra 
in the rest-frame wavelength range spanning [Oil] to Ha. A 
detailed outline of the spectral reduction pipeline is given 
in CoUess et al. (2001) and F99 which we follow with only 
minor modifications. 

The data that we have available for this analysis include 
all observations up until January 2001. At this stage a total 
of 111404 and 60062 galaxy spectra had been obtained in 
the SGP and NGP regions respectively. From these objects 
we remove those which did not receive accurate redshifts 
(Q < 3, CoUess et al. 2001), repeated observations and spec- 
tra with particularly low signal-to-noise ratio (< 10). This 
leaves us with 78994 (SGP) and 44937 (NGP) galaxies for 
use in our spectral analysis. Because of interference from sky 
emission and atmospheric absorption we have made a fur- 
ther restriction to our redshift range such that z < 0.15. This 
cut ensures that the Hq line in our spectra is not corrupted. 
Note that this is only a problem for the small redshift range 
0.15 < 2: < 0.17, however, we have also excluded the galax- 
ies with 2 > 0.17 in order to simplify our analysis. Imposing 
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Figure 1. The average spectrum of the M;, — 51ogj^Q(fc) > — 18 volume-Umited sample of galaxies. The main spectral features present 
are labelled. 



this cut and removing spectra observed in poor quality fields 
(Section 4.1) leaves us with 43449 (SGP) and 32140 (NGP) 
galaxies. A total of 75589 unique galaxies will therefore be 
used in our subsequent analysis and measurement of the LF. 



3 SPECTRAL TYPES 

The spectral classification presented here is based upon a 
Principal Component Analysis (PCA) of the galaxy spectra. 
PCA is a statistical technique which has been used with con- 
siderable success by multiple authors in the past (e.g. Con- 
nolly et al. 1995; Folkes, Lahav & Maddox 1996; Glazebrook, 
Offer & Deeley 1998; Bromley et al. 1998; F99) to deal with 
large multi-dimensional data sets. A detailed mathematical 
formulation of the PCA adopted here in is given in F99. Note 
that the most significant difference between this formulation 
and that used by other authors (e.g. Connolly 1995) is that 
our spectra have been mean-subtracted (Fig. hi) before con- 
structing the covariance matrix. This makes no substantial 
difference to the analysis since using other methods simply 
yields the mean spectrum as the first component. 

We note that throughout the remainder of this pa- 
per we will denote the eigenvectors (herein eigenspectra) as 
PCi,PC2 etc. and the projections onto these axes by pci, 
pc2 etc. 



3.1 Applying PCA 

PCA is a useful technique in that it allows us to easily vi- 
sualise a multi-dimensional population in terms of just a 



Table 1. The relative importance (measured as variance) of the 
first 8 principal components. 

Component Variance (%) Component Variance (%) 
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handful of significant components. It does this by identi- 
fying the components of the data (in this case the galaxy 
spectra) which are the most discriminatory between each 
galaxy. The significance of each component is measured in 
terms of its contribution to the variance over the sample and 
is determined in the PCA. This allows us to identify just the 
most significant components for future use. It is clear from 
such a formalism that any clustering in a space defined by 
the PCA is indicative of distinct sub-populations within the 
sample 

In terms of reduced dimensionality we find, after apply- 
ing the PCA to our galaxy spectra, that rather than using 
the original 738 spectral channels to describe each spectrum 
we can use only two projections and still retain two thirds 
of the total variance within the population. Indeed we can 
retain 50% of the total variance just from the first compo- 
nent (the total variance of the population includes noise and 
hence this is an underestimate of the true fraction of vari- 
ance in these projections). The significance of each succes- 
sive principal component drops off very sharply so these first 
few components are by far the most important (Table hi). 

Admittedly it is not clear a priori whether retaining 
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maximum variance over the spectral population also retains 
a corresponding maximum of physical information regard- 
ing each galaxy. However, in this paper, and also in previ- 
ous independent analyses, it is found empirically that the 
components dominating the spectral distribution do in fact 
relate well to physical attributes For example, morphology 
and other well known diagnostics such as the continuum 
slope (colour), the average emission/absorption line strength 
and the strength of the Ha line (see e.g. Ronen, Aragon- 
Salamanca & Lahav 1999). 

It is worth noting here that there are several limitations 
to using PCA which the reader should bear in mind. The 
most important of these is that it imposes linearity on the 
data set. This is certainly not a satisfactory approximation 
in all cases, however the general ability of the analysis to suc- 
cessfully reconstruct galaxy spectra with only the first few 
projections suggests that such an approximation is valid. We 
also note that in this analysis (as in previous work) we have 
restricted ourselves to deriving only the most general classi- 
fication scheme. In so doing we ignore several distinct sub- 
populations in our data set such as star-bursts and AGNs. 
In fact it should be possible to take another approach and 
to target the analysis on finding these sub-populations. 



3.2 Volume-limited PCA 

Due to the large number of galaxy spectra acquired so far in 
the 2dFGRS it has now become possible to consider volume- 
limited subsets of the galaxies for use in the PCA. This 
approach is much more optimal compared to considering the 
entire flux-limited data set in which the spectral population 
can be biased by relationships between the luminosity and 
type of galaxies. 

We find that using an absolute Mtj — 5\ogiQ{h) magni- 
tude limit of —18.0 (corresponding to a redshift of z ~ 0.1) 
gives a representative subset of the local population whilst 
still retaining 35,289 galaxies in the analysis, a very large 
fraction of the total. We have therefore chosen to perform 
the PCA on this subset of the population in order to deter- 
mine the principal components. It is then a trivial matter to 
use these components to derive the individual projections for 
each of the remaining galaxies outside of our volume limited 
sample. The average spectrum of this volume- limited sample 
is shown in Fig. fjl and the projections in the plane defined 
by the first two principal components are shown in Fig. 0. 



3.3 Interpreting the PCA 

We can identify general trends in the distribution of projec- 
tions shown in Fig. by considering the physical significance 
of the first two eigenspectra. We find that whilst the first 
eigenspectrum contains an approximately equal contribution 
from both the continuum and the emission/absorption line 
strengths, the second is dominated by the latter. Therefore 
by taking certain linear combinations of the two we can ei- 
ther maximise or minimise the contribution from the line 
features in the spectrum (Fig. H). Such a transformation of 
axes allows us to easily interpret the distribution of Fig. g 
as both a sequence in colour from red (bottom left) to blue 
(top right) and an orthogonal sequence from absorption (top 
left) to emission (bottom right). 
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Figure 2. Projections (pci and pc2) in the space defined by 
the first two eigenspectra derived from the PCA for the Mf,j — 
51og]^Q(/i) < — 18 sample of galaxy spectra. The trends illustrated 
have been derived from the eigenspectra by identifying the linear 
combinations that either maximise the infiuence of the absorp- 
tion/emission features or minimise them. 



3.4 Instrumental effects 

The 2dF instrument was designed to measure large numbers 
of redshifts in as short an observing time as possible. In so 
doing it has made very large projects such as the 2dFGRS 
possible. However, in order to optimise the number of red- 
shifts that can be measured in a given period of time, com- 
promises have had to be made with respect to the spectral 
quality of the observations. Therefore if one wishes to char- 
acterise the observed galaxy population in terms of their 
spectral properties care must be taken in order to ensure 
that these properties are robust to the instrumental uncer- 
tainties (see e.g. Lewis et al., 2001, for a more detailed dis- 
cussion) . 

The 2dF instrument makes use of up to 400 optical fi- 
bres with a diameter of 140^m (corresponding to 2.0 — 2.16" 
on the sky, depending on plate position). The quality and 
representativeness of the observed spectra can be compro- 
mised in many ways, the most significant of which are: 

• Astrometric errors in determining the actual position of 
the galaxy on the sky can result in the fibre being placed in 
a position that is not optimal for light collection. In general 
the position of a galaxy centre is known to within 0.4" . 

• Given a known galaxy location on the sky the robotic 
positioner will place the fibre to within 20/jm. The RMS 
error in placing this fibre is found to be ll^m (0.16"). 

• The observed colour (continuum slope) of the galaxy 
depends on the positioning of the fibre aperture. This is due 
to a feature of the 2dF corrector lens design which results in 
relative spatial displacements between different components 
of a galaxy's spectrum. These displacements can be as large 
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Figure 3. The first two eigenspectra are shown in the top two panels and the linear combinations which either minimise (a) or maximise 
(b) the effect of emission/absorption features are shown in the bottom two. 



as 1" (at the 20' radius) between the incident blue and red 
components of the observed spectrum. This together with 
the previous two uncertainties in placing the fibre aperture 
will result in uncertainties in the continuum measurement 
which make the spectra very difficult to accurately flux cal- 
ibrate. 

• The apparent size of a galaxy on the sky can be larger 
than the fibre aperture. The positioning of this aperture 
can therefore result in an unrepresentative spectrum due to 
spectral gradients within the galaxy. 

• Atmospheric dispersion (due to imperfect correction) 
and the differential atmospheric refraction during an expo- 
sure each result in errors of the order of < 0.3". 



tained in the emission/absorption line strengths should be 
relatively robust. 

Perhaps our greatest concern regarding the represen- 
tativeness of the spectra is due to the limited size of the 
aperture with respect to the galaxy. Surprisingly, we have 
found it difficult to isolate any systematic bias from this ef- 
fect. Of course this effect will be substantially diluted in the 
observed spectra due to the significant seeing present at the 
Anglo-Australian Telescope in Siding Springs, which is of 
the order of 1.5" — 1.8". In addition the presence of differ- 
ential atmospheric refraction may also assist us. We discuss 
the impact of this 'fibre-aperture bias' on our results in more 
detail in Section 5.3. 



The source of the chromatic dispersion mentioned above 
is due to a design feature of the 2dF corrector lens which 
gives us our extended (2 degree) field of view. This disper- 
sion, combined with the random positioning errors, results 
in uncertainties in calibrating the continuum of each galaxy 
spectrum. For this reason we will not make use of the in- 
formation contained in the continuum of each galaxy in our 
subsequent analysis. However, we will make use of the con- 
tinuum calculated by averaging over a sufficiently large en- 
semble of galaxy spectra since this will be much more robust 
(see Section 4.2 for further details). 

The random errors in positioning the fibre aperture 
could also result in uncertainties in measuring the line- 
strengths of emission features such as Ha since we would 
expect the strength of this emission to vary spatially across 
the galaxy. However, because the random positioning errors 
are much smaller than the scale lengths of the galaxies this 
isn't a large problem. For this reason the information con- 



3.5 Tj parameterisation 

We are presented with several options in attempting to de- 
rive a classification for the observed 2dF galaxies based upon 
their spectra. However, there are several issues that one must 
first consider: 

• The distribution of galaxies in the {pci,pc2) plane is 
smooth (although slightly bimodal), representing a contin- 
uous sequence from absorption to emission and from red to 
blue. There is no evidence to suggest any division into dis- 
tinct classes of galaxy spectra. 

• The instrumental errors result in unstable continuum 
measurements. On the other hand small scale features such 
as line-strengths are relatively unaffected. 

• The relationship between morphology and spectral 
properties is not very well understood and hence attempting 
to anchor the classification using a morphological training 
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Figure 4. The observed distribution of spectral type measured by 
Tj. Also shown are the four divisions we use to divide the datasct 
(see Fig. M for mean spectra) . The bottom panel shows the corre- 
lation between r) and morphological type using a training set of 
galaxies taken from the Kennicutt Atlas (Kennicutt 1992). 



Figure 5. Dispersion in the PCA projections (pci and pc2) and 
in the spectral classification rj. These distributions are well fit by 
Gaussians (shown here) and were calculated by taking the differ- 
ences in these quantities between ~ 2000 repeated spectral obser- 
vations. Also shown is the dispersion in the orthogonal projection 
to Tj, which quantifies the continuum slope in the spectra. 



set (e.g. Connolly et al. 1995; F99) would be premature. A 
more robust and quantifiable measure of the galaxy type is 
required. 

By projecting the pci and pc2 components of each galaxy 
onto the linear combination which maximises the effect of 
emission/absorption line features we are in effect high-pass 
filtering the spectra. Thus we would expect (and indeed we 
find) this projection to be relatively stable to uncertainties 
in continuum measurements. 

By using this projection we are determining a measure 
of the average emission/absorption line strength of a galaxy 
which is easily quantifiable and robust. In addition this pro- 
jection is also representative of the spectral sequence of the 
galaxy population since it is composed of the two most sig- 
nificant principal components (representing 66% of the total 
variance over the population). 

We therefore choose to adopt this projection, which we 
shall denote r), as our continuous measure of spectral type 



rj = a pci — pc2 



(1) 



The value of a which maximises the emission/absorption 
features is essentially identical to that we would find from 
identifying the most stable projection in the (pci,pc2) plane 
between repeated pairs of spectral observations. Using this 
method we find a — 0.5 ± 0.1. 

In Fig. the distribution of the rj projections is shown 
for the galaxies observed to date in the 2dFGRS. Also shown 
in the same figure is the ry-morphology relation for a sam- 
ple of galaxies from the Kennicutt Atlas (Kennicutt 1992) . 
Comparing the two data sets shows that there is a corre- 
spondence between the sequence of r/ and that of morphol- 



ogy. Note that this correspondence can only be treated as 
an approximation since our sample of Kennicutt galaxies is 
not complete and only represents a few high signal-to-noise 
ratio spectra with well determined morphologies. However, 
the trend is clear. 

It is important to have an understanding of the un- 
certainties in the principal component analysis as this will 
affect the reliability of our classification. Using a sample of 
(~ 2000) repeated spectral observations we can attempt to 
quantify these uncertainties. The differences in the PCA pro- 
jections (pci and PC2) and in the spectral classification be- 
tween repeated pairs is shown in Fig. H. Note that for clarity 
we show here the Gaussian fits made to these error distri- 
butions. It is clear from this figure that the uncertainty in 
the pci and pc2 projections is substantial (la dispersions of 
2.9 and 1.7 respectively), however, this dispersion is greatly 
reduced by using the rj parameterisation (Icr = 0.7). 

In addition to correlating with morphology, rj has also 
been shown to be strongly correlated with the equiva- 
lent width of Ha in emission line galaxies. Fig. (Bland- 
Hawthorn et al., in preparation). Hence rj can be interpreted 
as a measure of the current star-formation present in each 
galaxy. 

Note that having a continuous measure of spectral type 
allows us to take greater advantage of having such a large 
data set in which much of the detail in the spectral sequence 
would be smeared out by rough divisions (as in the case 
of morphological segregation). However, for the purposes of 
this preliminary analysis we choose to split the spectral se- 
quence into four broad bins ('types') as shown in Fig. W. 
The divisions we have made are not determined by the mor- 
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Figure 6. The EW(Ha) is very tightly correlated to rj for emis- 
sion line galaxies (see Bland-Hawthorn et al., in preparation, for 
further details). 



the STY LF estimator (Section 4.3). In the case of appli- 
cations where the redshift limits are explicitly set this com- 
pleteness ratio is adjusted to take account of this, although 
this is only a small correction. We find that the classification 
completeness can be well fit by a function of the form. 



C(b, 



- exp[/3(6,, - 7)] 



(3) 



We determine the best fit parameters by weighting each 
bin according to the number of galaxies that contributed 
to it and performing a non-linear least squares fit, yield- 
ing (q,/3,7) = (0.96,1.51,20.73). A weight for each galaxy can 
then be determined from the inverse of this quantity. 

In addition we note that this parameterisation can be 
significantly skewed in observed fields with low overall com- 
pleteness. So before calculating these parameters we have 
imposed a minimum completeness threshold for each ob- 
served field of 90%. This additional cut reduces the number 
of galaxies to be used in our calculations by approximately 
8,000 (as described in Section 2). 

The completeness of our sample may also be infiuenced 
by other selection biases which vary systematically with red- 
shift or spectral type e.g. surface brightness. Of course it is 
not possible to determine a simple correction for these ef- 
fects since we do not know the redshift or spectral type of 
a galaxy until it has been successfully observed. We expect 
these additional effects to be small, they will be discussed 
in more detail in a forthcoming paper (Cross et al., 2001). 



phological training set (as in previous work) but rather by 
the shape of the 77-distribution itself. The peak in the distri- 
bution arises from the degeneracy in spectral line-features 
between elliptical and early-type spiral galaxies. We sepa- 
rate this peak from the shoulder, which we then divide into 
two. Another cut is then made to separate the tail of the 
distribution which will be dominated by particularly active 
galaxies such as starbursts and AGN. The average spectrum 
of each type is shown in Fig. H. 



4 LUMINOSITY FUNCTIONS 

4.1 Completeness 

We find that at a depth of 6j = 19.2 the density of galax- 
ies is 149.4/0° over our targeted area of the sky (Norberg 
et al., 2001). In order to avoid fibre collisions it is not pos- 
sible to assign a fibre to every galaxy in a given field. We 
take into account this configuration completeness (relative 
to the parent catalogue) of 93% to calculate the effective 
area of the observed galaxies, thereby providing the overall 
normalisation for our LF estimates. 

In addition to the configuration completeness it is also 
necessary to correct for the classification incompleteness: we 
chose to exclude spectra with Q < 3 and signal-to-noise ratio 
< 10, which do not yield an accurate value of rj. We assume 
that this completeness can be parameterised only in terms 
of the apparent 6,, magnitude and define 



C(b,, 



^Galaxies that could be classified 



(2) 



^Galaxies targetted 

This ratio is independent of redshift and is appropriate for 
applications based on conditional probabilities p[M\z) e.g. 



4.2 iv-Corrections 

Here we develop a self-consistent method for i('-corrections 
based on the observed 2dF spectra. This is a significant 
improvement compared with previous analyses where the 
_K'-corrections were anchored to standard SEDs (e.g. Pence 
1976) through the relation between morphology and spectral 
type. 

The main obstacle in calculating A'-corrections directly 
from the observed spectra is the large uncertainty associated 
with measuring the continuum (Section 3.4). If this were 
not an issue then the calculation of the ii'-corrections would 
be straightforward, since both the rest-frame and observed 
wavelength range of the 2dFGRS spectra span the range of 
the foj passband where the transmission is above 5% of the 
passband peak (3710 - 5520A). 

On a galaxy-by-galaxy basis the continuum uncertain- 
ties caused by the chromatic distortion in the telescope op- 
tics coupled with the small random fibre positioning errors 
will induce random errors in the _K'-corrections. We are as- 
sisted by the fact that there are only three substantial sys- 
tematic effects which might affect our calculations: 

(i) Measuring the &j passband. (Hewett & Warren, pri- 
vate communication). 

(ii) Determining the 2dF system response (Lewis et al., 
2001). 

(iii) Aperture bias (Sections 3.4 and 5.3) due to the rel- 
atively small size of our fibre aperture compared with the 
spatial extent of low-redshift galaxies. 

Of these three the 2dF system response is by far the least 
well determined. The system response is calculated using 
photometric standard star observations with the 2dF instru- 
ment setup and accounts for the relative photon capturing 
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Figure 7. The average 2dFGRS spectrum of each spectral type is shown. Each spectrum is in units of counts per pixel and has been 
normalised to have average counts of unity. 



efficiency of the instrument over different wavelength bands 
(corrected for seeing). The exact nature of this response is 
still in the process of being fully determined. For the mean- 
time we adopt a small correction to the previously measured 
standard-star system response in order to correct the aver- 
age galaxy survey observation, 1997-2001, (Baldry et al., 
2001). 

The aperture bias mentioned in the third point can have 
quite a significant impact on the observed 2dFGRS spectra 
since for low-redshift galaxies we will systematically sam- 
ple more of the (redder) inner bulge light in our spectra. 
However, this will not lead to significant errors in the ab- 
solute luminosities since at low redshifts, where the aper- 
ture bias can be substantial, the A'-corrections are small 
in absolute terms; conversely, at higher redshifts where the 
.ft'-corrections are large, the aperture bias is minimal. We 
have conducted several tests to estimate the significance of 
aperture bias on our spectra and we have found it to be 
negligible beyond z ~ 0.1, particularly due to the relatively 
poor seeing at the AAT which effectively 'smears' out spec- 
tral gradients in each galaxy. In addition we find that it will 
only have a significant impact for particularly low-redshift 
galaxies. 

Each individual galaxy spectrum has large uncertainties 
associated with its continuum, so to calculate meaningful 
Tf- correct ions we adopt here a practical approach. Having 
divided the galaxies into their respective classes we proceed 
to calculate the average spectrum for each class. Since each 



class contains many thousands of galaxies it is reasonable 
to expect any unphysical features present in the spectra to 
average out (assuming we have accounted for the above sys- 
tematic uncertainties). We find from repeated spectral ob- 
servations that it is equally likely for a continuum to be 
bent one way as the other, so assuming we have averaged 
over a sufficient number of galaxy spectra it is reasonable to 
assume that these effects will balance out. Therefore each 
average spectrum can be considered to be a characteristic 
galaxy template that reflects the true average SED of that 
ensemble of galaxies. We can now proceed to calculate K- 
corrections for each of these four templates. 

Note that we apply these if-corrections to each galaxy 
individually to calculate its absolute magnitude. However, 
on a galaxy-by-galaxy basis these derived absolute magni- 
tudes will not be optimal. Rather, they should be inter- 
preted as part of the ensemble to which they belong. We 
should therefore only apply these /('-corrections to appli- 
cations which involve comparing ensemble properties (e.g. 
LFs) or alternatively to calculating a property of the entire 
data set. 

We have fit quadratic forms to the 6j A'-corrections 
given by the following expressions; 



fci(2) = 2.62 + 4.32 
k2(z) = 1.9z + 2.2z^ 
kz{z) = 1.3z + 2.0z'^ 



(4) 
(5) 
(6) 
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Table 2. The mean shifts between the 2dFGRS 6j magnitudes 
and the CCD magnitudes per spectral type. The stated uncer- 
tainties are la. 



Type Mean shift (mag) Uncertainty 



-0.009 

-0.033 

0.024 

0.09 



0.005 
0.004 
0.006 
0.01 



ki{z) =0.92 + 2.3z 



(7) 



Where ki refers to the iC'-correction of type i and is valid over 
the redshift range < z < 0.2. The average iC'-correction 
over all spectrally classified galaxies is given by fcav(z) = 
l.'dz + 2.7z^ . This expression can be used in lieu of the above 
or instead could be applied to non-classified galaxies which 
require ii'-corrections. 

The most substantial source of errors in our K- 
corrections will arise from the uncertainty in the 2dF system 
response. Taking this into account we estimate the uncer- 
tainty in our ii'-corrections to be of the order of 20%. 

Norberg et al. (2001) have also calculated is'-corrections 
for the 2dFGRS, independently of the observed galaxy spec- 
tra. They find very similar 7^-corrections, which are consis- 
tent with our own to within the stated uncertainty. 



Figure 8. The shift between the 2dFGRS 6j magnitude and that 
derived from a CCD photometric survey (Cross et al., in prepa- 
ration) is shown for ~ 3000 galaxies. It can be seen that the shift 
varies systematically with surface brightness and so care must be 
taken when dividing the sample. 



4.3 Surface brightness effects 

The total 6,, magnitudes used in this analysis have been de- 
rived from plate scans (Maddox et al. 1990) and so they will 
be more susceptible to certain systematic effects than those 
derived from CCD data. In order to improve the accuracy of 
our magnitudes several re-calibrations of the 6j magnitudes 
have been made using CCD photometry from overlapping 
fields (see Norberg et al., 2001). 

Unfortunately some small offsets still remain. One such 
effect which may be important to the analysis of LFs per 
spectral type is an apparent shift in 6j magnitude with sur- 
face brightness (Cross 2001; Blanton et al., private commu- 
nication) . Fig. y shows the offsets between the 2dFGRS 6j 
magnitudes and those derived from a CCD photometric sur- 
vey (Cross 2001; Cross et al., in preparation) versus surface 
brightness, a trend is clear in that high surface brightness ob- 
jects tend to be systematically fainter in 2dFGRS 6j. These 
offsets are most likely due to a saturation effect in the (plate 
derived) 2dFGRS magnitudes. 

These small offsets will have a negligible effect on the 
calculation of the overall LF but may be significant when 
one divides the 2dFGRS sample in a way which is related 
to surface brightness. Fig. shows how these shifts depend 
upon our spectral classification, it can be seen that the re- 
lationship is surprisingly weak. 

The mean shifts we calculate per type are however 
significant compared to their estimated uncertainties (Ta- 
ble. 0), particularly for the most late- type galaxies. We cor- 
rect the galaxy magnitudes for each spectral type in order 
to account for these shifts in the analysis that follows. 
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Figure 9. The bj magnitude error versus spectral type (tj) is 
shown. It can be seen that there is a weak correlation between 
the two quantities. 



© 2001 RAS, MNRAS 000, 0-0 



10 Darren S. Madgwick et al. 



4.4 Estimators 

We assume a flat homogeneous Universe witii a uniform 
Hubble flow {h = iifo/lOOkms"^ Mpc~^), a cosmological 
constant JIa = 0.7 and mass density parameter f2ni = 0.3. 
In computing tiie LFs we use botli the step-wise maximum 
likelihood (SWML) method (Efstathiou, Elhs & Peterson 
1988, hereafter EEP88) to make a non-parametric fit to the 
LF and the STY method (Sandage, Tammann & Yahil 1979) 
to calculate the maximum-likelihood Schechter function fit 
to the LF, 



<f>{L)dL - 



<h) 



exp 



L_ 
I* 



dL 
U 



(8) 



Both of these estimators are discussed extensively in the lit- 
erature (e.g. EEP88; Willmer 1997). In addition we take into 
account the presence of Malmquist bias in our STY fit by 
assuming that the errors in the 6j magnitudes are Gaussian 
{aM = 0.15, see CoUess et al. 2001). The observed LF wUl 
then be a convolution of the true LF with the distribution 
of magnitude errors, 



<^obs(M) 



2ticjm 



<j,(M')e-^^'-^'^"^'^MdM' . (9) 



Neither method predicts the normalisation of the LF so 
we estimate this independently by fitting to the projected 
counts, 149.4/D° (Section 4.1). In so doing we make the as- 
sumption that the type fractions in the parent catalogue are 
identical to those of the classified galaxies (over the redshift 
interval 0.01 < z < 0.15). 



5 DISCUSSION 

5.1 Results 

Having performed the calculations outlined above we find 
that the Schechter function provides an inadequate fit of the 
LF calculated over the entire range in Mtj magnitude, espe- 
cially for our most passive and active star-forming galaxies 
(Types 1 and 4 respectively). In order to determine a better 
fit to the LF of the passive (Type 1) galaxies we have applied 
the following simple procedure: First we fit a Schechter func- 
tion to the LF of the bright galajcies (Mj,, — 51ogj^Q(/i) < —17) 
using the STY method. This provides an adequate represen- 
tation of the data over the brightest ranges in magnitude. 
Having calculated this fit we then found that the residual 
LF at the faint end could be well fit by a function of the 
form 



0ros(M6,) = lo["+''*^''.i', 



(10) 



particularly for the magnitude range —17 < Mbj — 
51og]^Q(/i) < —14.5 where the error-bars are small. We calcu- 
lated a and b using a least-squares fit to the residual created 
by subtracting the STY fit from the SWML points at the 
faint end. The LF at the faint end can then be described by 
the combination 



(j){M) = (l>sch{M; a, M' ,<!>*) + (j),^4M; a,b) , 



(11) 



where (a, b) = (4.7, 0.50). The uncertainties in these param- 
eters are quite large (±0.9, ±0.06), an obvious consequence 
of the large errors in the SWML points at these magnitudes. 



We estimate the scatter in our derived Schechter param- 
eters by performing a large number of random realisations 
drawn from the estimated LFs. In addition to this sampling 
variance we also take into account the uncertainties in our 
_R'- corrections (20%) and in the number density of the galax- 
ies (7%). The quoted errors are at the la level. 

The two most notable features of our derived LFs are 
the systematic variation of the faint end slope (a) and the 
characteristic magnitude (M*) with spectral type. The faint 
end slope steepens significantly from a = —0.54 for the most 
passive star- forming galaxies (Type 1) to a = —1.50 for 
the most active galaxies (Type 4). On the other hand the 
characteristic magnitude becomes systematically fainter as 



we go from the most passive {A4* — 51ogj^Q(/i) 



-19.58) 



to active {Al* — Slog^oCi) = —19.15) star-forming galaxies. 
These results agree with previous measurements in that the 
late- type galaxies are systematically fainter than their early- 
type counterparts. 

Note that our results do not change much if we exclude 
our correction for the errors in our 6j magnitudes (AA/* = 
0.05 - 0.06, Aq = 0.02 - 0.05). 

The fact that the Schechter function is not a good fit to 
the data over the entire M^j magnitude range is mostly due 
to an over- (Type 1) or under- (Type 4) abundance of faint 
objects relative to bright objects. This is particularly true 
of the passive (Type 1) sample of galaxies for which there is 
a very significant increase in the predicted number density 
of faint galaxies. This feature has been present in previous 
analyses (e.g. F99) however the small sample size has meant 
that only a statistically insignificant number of galaxies have 
contributed. Hence previous studies could not determine if 
this feature was real or a consequence of the small volume 
being sampled at these magnitudes. For the flrst time here 
we show signiflcant evidence for the presence of a substantial 
dwarf (passive star-forming) population with 142 galaxies 
having Altj — 51ogiQ(/i) magnitude fainter than —16.0 {z = 
0.027). Of course these objects will tend to have low signal- 
to-noise ratio spectra so some degree of contamination in 
this sample should be present, particularly from erroneous 
redshift determinations. However, these preliminary results 
look promising and we expect more detailed analyses to be 
forthcoming in the near future. 



5.2 Previous results 

F99 has performed similar calculations of the LF per spectral 
type using a subset of the data presented here. In Fig. hll we 
compare their results with those presented here. Note that in 
order to make this comparison we have recalculated the LFs 
assuming (for consistency) an Einstein de-Sitter Universe 
(fi = 1, A = 0), and we have reduced the 5 classes of F99 
to our 4 types by taking linear combinations of their LFs 
which give the same fractions of galaxies per type as our 
own. It can be seen from the flgure that our results are in 
good agreement with those of F99 except for some small 
systematic shifts. These differences are due to four factors 
which we believe have significantly improved the accuracy 
of our analysis. 

(i) We have used a considerably larger sample of galaxies 
in our calculations. 
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Figure 10. The luminosity functions derived per spectral type. The top left box shows the LF derived from all the classified galaxies, 
and the top right box shows how this (solid line) relates to the summation of the fitted LFs per type (dashed line). For each of the 
remaining four panels we show the SWML points and Schechter fits to the LF per spectral type. The Schechter function has been fit 
over the entire range in M{,j magnitude for Types 2, 3 and 4, and only out to Mn^ — 5\ogig{h) = — 17 for Type 1. For Type 1 we show 
the best fitting Schechter function for the bright galaxies (dotted line) and the combination of this with our analytic form for the faint 
end residuals (solid line). Our assumed cosmology here is f2m=0.3, Oa=0.7. 



(ii) We have adopted self consistent iC-corrections, result- 
ing in a small change in the predicted M* and (j>* ■ 

(iii) Taking into account recent advances in our CCD cal- 
ibrations we have been able to update the catalogue 6j mag- 
nitudes (Colless et al. 2001; Norberg et al., 2001). This has 
also resulted in a small change in the number counts and 
hence 0*. 

(iv) We have corrected for saturation effects in high sur- 



face brightness objects (Section 4.3). This has a particularly 
large impact on our most late-type (Type 4) galaxies. 

The most similar independent analysis of the LF per 
spectral type has been conducted by Bromley et al. (1998), 
using the Las Campanas Redshift Survey. We also compare 
these results with our own in Fig. hU Again, in order to com- 
pare the 6 clans with our 4 types we have had to combine 
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Table 3. Schechter function parameters derived in our analysis per spectral type using galaxies with Mj < Mn^. The 



errors given are itlcr. The normalisation, <j>* , of the Schechter function is in units of 10 
is the luminosity density in units of IO^Lq/iMpc"^ (^^b^Q ~ 5-30). 



'ft^Mpc-'^ and PL = 0*L*r(Q + 2) 



Sample 


■n 


Mim 


- 5 logir 


(h) 


# Galaxies 


Af*-51ogioW 


a 


r 


PL 


All 


- 




-13.0 




75589 


-19.79 ±0.04 


-1.19 ±0.01 


15.9 ±1.0 


19.9 ±1.4 


1 


V< -1.4 




-17.0 




27540 


-19.58 ±0.05 


-0.54 ±0.02 


9.9 ±0.5 


7.8 ±0.5 


2 


-1.4<r;< 1.1 




-13.0 




24256 


-19.53 ±0.03 


-0.99 ±0.01 


7.2 ±0.4 


6.1 ±0.4 


3 


1.1 <r; <3.5 




-13.0 




15016 


-19.17 ±0.04 


-1.24 ±0.02 


5.0 ±0.3 


3.7 ±0.3 


4 


r;>3.5 




-13.0 




8386 


-19.15 ±0.05 


-1.50 ±0.03 


2.4 ±0.2 


2.6 ±0.3 



their LF determinations in a way which yields the same frac- 
tions of galaxies per type as our own. In addition we have 
applied the transformation M^^ = M,- + 1.1 (Lin et al. 1996) 
to convert the LCRS Gunn-r Magnitudes to fej. 

The trend of the galaxy LFs to become progressively 
fainter as one moves from passive to active star-forming 
galaxies is clear for both sets of results. However there ap- 
pears to be some significant disagreement particularly in the 
faint end slope between the two LF determinations. Given 
the poor fit of the Schechter function to the faint end of the 
LF we would not expect these two sets of results to necessar- 
ily agree within their stated errors. The large differences ob- 
served are most probably due to the different selection effects 
(e.g. surface brightness and foj versus r selection) which can 
be particularly prominent for fainter galaxies, hence adding 
more uncertainty to a than can be estimated from a stan- 
dard error analysis. 

Our total luminosity density (again assuming an Ein- 
stein de-Sitter Universe) is within 4% of the figure obtained 
in F99 and is 30% lower than the figure quoted by Blanton 
et al. (2001) for the Sloan Digital Sky Survey early data 
release. The reasons for this discrepancy are discussed in 
Norberg et al. (2001), but in summary we believe that the 
Blanton et al. figure is too high, partly owing to large-scale 
structure in their small initial sample, and also through the 
use of an inappropriate colour equation to predict 2dFGRS 
magnitudes. 



5.3 Aperture effects and evolution 

The spectra observed as part of the 2dFGRS are acquired 
through optical fibres with a diameter of 140/im, correspond- 
ing to between 2" and 2.16" on the sky (depending on plate- 
position, see Lewis et al., 2001). In our assumed cosmology 
(ilm=0.3, riA=0.7) this translates to an aperture width of 
approximately 2.8h~^ kpc at the survey median redshift of 
z — 0.1. The fact that fitted galaxy brightness profiles have 
scale lengths of the order of 3-4 kpc (Binney & Merrifield 
1998) raises the question of how representative the observed 
spectra are of the galaxies as a whole. However, with typ- 
ical seeing (1.5" — 1.8") and random fibre positioning er- 
rors (RMS 0.3") the 2dF fibre-aperture will in fact sample 
a much larger portion of each galaxy. 

If the effective fibre apertures were significantly smaller 
than the galaxies over a range of redshifts then the varia- 
tion in this aperture with redshift would result in system- 
atic biases in the determination of the LFs. These so-called 
'aperture effects' are difficult to quantify without an in depth 




Figure 11. Comparison with the type-dependent LFs calculated 
by Folkes et al. (dashed lines) and Bromley et al. (dotted lines). 
The solid lines show our LFs assuming (for consistency) an Ein- 
stein de-Sitter Universe. Note that the Type 1 LF has only been 
fit to Mi,j — 5logiQ{h) < —17. The disagreement between the 
different sets of determinations may be due to the poor fit of 
the Schechter function to the LF or to different selection effects 
between the surveys. 



knowledge of the spatially resolved spectral properties of the 
targeted galaxies. 

Kochanek, Pahre & Falco (2000) recently suggested a 
test to estimate the infiuence of this fibre-aperture bias on 
the determination of type dependent LFs (j>i{M). The test 
involves comparing the observed fractional number counts of 
galaxies in each class to the fractional number counts pre- 
dicted by our derived LFs. Kochanek, Pahre & Falco found 
that for redshift surveys classified morphologically the two 
were generally consistent but that for surveys classified spec- 
trally using fibre-aperture spectra the two were very incon- 
sistent. We have performed this test on our data and found 
no substantial inconsistencies (Fig. 113). With the exception 
of small systematic offsets the observed type fractions were 
found to agree with those predicted by the LF to within 2cr 
over the entire range in 6j magnitude. We note that as the 
number counts predicted by the LF are normalised by con- 
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Figure 12. A test recently proposed by Kochanek, Pahre &; Falco 
(2000). We show observed (squares) and predicted (lines) type 
fractions for the 2dF galaxies used in our calculations, where the 
predicted fractions have been derived from integrating our LFs. 
The two are consistent within the stated (Poisson) errors. 



Figure 13. Comparison between the observed N{z) distributions 
(histograms) and the predicted N{z) for each spectral type (6j < 
19.3). Note that the A'^(2) distributions of the first two types are 
recovered very well but that significant deviations exist for the 
last two (most active star- forming) types. 



struction to the observed number counts this test is only nec- 
essary and not sufficient. For this reason more sophisticated 
tests such as comparisons with detailed models are neces- 
sary in order to isolate and analyse the competing effects on 
the data (e.g. positional errors, seeing and evolution). 

We show the observed and predicted N{z) distributions 
for each spectral type out to bj < 19.3 in Fig. uA It can be 
seen that the observed N{z) is well recovered by the LF 
for the first two types. There appears to be a significant 
over abundance of high-2 galaxies for the two highest star- 
forming types as well as a slight under abundance of low-z 
galaxies. Similar deviations in abundance could be explained 
by aperture effects. However, if this were the only explana- 
tion then we would also observe a significant over abundance 
of \ow-z passive star-forming galaxies, which is not the case. 
A much more reasonable explanation for these deviations 
between the predicted and observed N{z) distributions is 
the presence of evolution (see e.g. Lin et al. 1999). Because 
these two effects affect the observed galaxy population in 
similar ways they will be difficult to dis-entangle in more 
detailed analyses. For the present we note that as far as the 
determination of the LF in broad bins of spectral type is 
concerned both of these biases have only a small effect on 
our calculations. 



6 CONCLUSIONS 

In this paper we have applied a method of spectral clas- 
sification, ri, based on the the relative emission/absorption 
line strength present in a galaxy's optical spectrum. This 
classification correlates well with morphology and we inter- 



pret it as a measure of the relative current star-formation 
present in each galaxy. Based on this classification we have 
divided the galaxies observed to date in the 2dFGRS into 
4 spectral types and calculated the LF for each. Our re- 
sults are consistent with previous determinations in that the 
faint end slope systematically steepens and the characteris- 
tic magnitude becomes fainter as we move from passive to 
active star-forming galaxies. In addition we have been able 
to show that the Schechter function is a poor fit to the LF 
over a large range of Mb^ magnitudes, perhaps as a result of 
sub-populations within the data. The calculations presented 
in this letter represent the most accurate determinations of 
the optical fej galaxy LF to date. 

As r; is a continuous variable the next step in quantify- 
ing the variation of the LF with spectral type is to develop 
a bivariate formalism 0(Af, j;). Such a result would be ex- 
tremely valuable in that it would quantify the variation in 
the Schechter parameters with galaxy type, a{ri) and M*(r]), 
allowing the results to be applied much more objectively to 
such things as semi-analytic models from which spectra as 
well as other galaxy properties can be derived (e.g. Slonim 
et al. 2001 and references therein). We can then also extend 
this approach of comparing the observed LFs with model 
LFs by mimicking the observational selection effects (e.g. 
aperture bias). 

One of the key advantages of our spectral type, -q, is that 
it can easily be applied to any other set of galaxy spectra 
which occupy the relevant wavelength range. This will be 
invaluable when comparing results from the many different 
galaxy surveys currently underway (e.g. SDSS, VIRMOS, 
6dF etc.). 
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The 2dF Galaxy Redshift Survey has made available to 
the community 100,000 galaxy spectra and redshifts through 
their website, http://www.mso.anu.edu.au/2dFGRS/. Ad- 
ditional information regarding the spectral classifica- 
tion and LFs presented here can be obtained from 



tittp : //www . ast . cam . ac. uk/^twodf grs/ 
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